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We analyse the possibility to create two-mode spin squeezed states of two separate spin ensembles 
by inverting the spins in one ensemble and allowing spin exchange between the ensembles via a near 
resonant cavity field. We investigate the dynamics of the system using a combination of numerical 
and analytic calculations, and we obtain squeezing for a wide range of parameters. We also investigate 
the transfer of the squeezing properties to the cavity field and to an output mode from the cavity. 
Finally, we investigate how the squeezing is affected by effects of inhomogeneities which would be 
present in solid state implementations of the spin ensembles. 
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I. INTRODUCTION 

Ensembles of identical quantum systems allow identifi- 
cation of collective degrees of freedom, which may couple 
strongly to, e.g., radiation fields and thus operate as 
efficient interfaces for preparation, manipulation and stor- 
age of quantum states of light pp. We will in this paper 
study the prospect of creating squeezed quantum states 
between two separate ensembles of effective two-level sys- 
tems which both couple to a single quantized field mode 
of optical or microwave radiation in a cavity. 

The basic idea of our proposal is to initialize the first 
ensemble with all the two-level systems in the ground 
state and the other ensemble with all systems in the ex- 
cited state. The coherent transfer of one excitation from 
one ensemble to the other, mediated by a cavity photon, 
conserves the energy and total number of excitations in 
the system, but if we redefine the labelling of the ground 
and excited states of the inverted ensemble, the excita- 
tion transfer is formally equivalent to a simultaneous pair 
excitation of the ensembles. For a high degree of po- 
larization, collective ensembles of two-level systems are 
effectively described as harmonic oscillators, where the 
number operator counts the number of excited two-level 
systems, while the position and momentum quadrature 
operators describe collective observables associated with 
the coherences in the two-level systems or, equivalently, 
the difference in occupation of superposition states of 
the system. In a spin language, the number operator is 
equivalent to the vertical component of the collective spin, 
while the quadrature operators measure the horizontal 
collective spin components. Correlated pairwise excita- 
tions of two harmonic oscillator modes are known in the 
non-degenerate optical parametric oscillator (OPO) in 
quantum optics, where it leads to squeezing and EPR 
entanglement of optical fields [2], and it is the purpose of 
this manuscript to investigate the accomplishments of the 
similar process in ensembles of two-level systems. 

We have two physical systems in mind, as illustrated 
in Fig.l. In part a) of the figure, we sketch a cavity 
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formed by a transmission waveguide for microwave fields, 
and we indicate the location of two solid state ensem- 
bles with a large number of spins interacting resonantly 
with the cavity field. Numerous experiments [3H7] have 
demonstrated the strong coupling between a cavity field 
and spin ensembles, consisting of different dopant ions 




FIG. 1. Two physical set-ups for collective spin squeezing 
of two- level ensembles, a) Two ensembles a and 6 of two- 
level systems are prepared in opposite eigenstates, illustrated 
with spin directions up and down. The spins are situated 
in a transmission waveguide cavity which can be tuned into 
resonance with the spin transition, b) Two ensembles of atoms 
or ions with pairs of ground states and optically excited states 
are coupled to a classical laser field with Rabi frequency Q 
and to a quantized optical cavity field with coupling strength 
X, such that Raman transitions between the ground states 
are accompanied by the creation or annihilation of a photon 
in the cavity mode. In both physical implementations the 
idea is to prepare each ensemble in one of the two internal 
states, enabling correlated transitions in the two ensembles, 
by exchange of cavity photons. 
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with both electronic and nuclear spin degrees of freedom, 
and NV centers in diamond. Following a recent proposal 
for hybrid quantum computing [8], transfer of quantum 
states was recently demonstrated between an ensemble 
of NV centers and a transmon qubit via the cavity field 
in a similar setup|9j. Our present work is related to, 
and supplements, recent theoretical proposals to generate 
squeezed [TUJ and entangled states [TTJ rjj] of NV-center 
electron spins. 

In part b) of the figure, we sketch an optical cavity 
in which a collection of atoms or ions are trapped and 
interact with the cavity field and with a classical control 
field. Two-photon Raman transitions effectively imple- 
ment atomic ground state changes associated with ab- 
sorption and emission of single cavity photons. Strong 
collective coupling to an optical cavity has been observed 
in recent experiments with ions [13 and neutral atoms 
[141117] . Spin squeezed and entangled states of atomic 
ensembles have applications in quantum metrology and 
quantum information protocols [18-21 ], and protocols for 
their generation via cavity mediated interactions, different 
from the one proposed here, have been proposed [2"2l |2"3"] . 

In this work the ensembles will experience a strong 
collectively enhanced coupling to the cavity and, as in [TUJ, 
the non-classical correlations are established by inverting 
the population in one of the ensembles such that all the 
spins are excited and spin flips occur accompanied by 
excitation of spins in the other ensemble. 

In section [TT] we describe the physical systems and the 
theoretical model that we will use in our calculations. In 
Sec. |III| we confirm the squeezing for both short and long 
times by a combination of numerical calculations and 
analytical arguments. In Sec. |IV| we consider the limit 
where the detuning of the cavity is large with respect to 
the degenerate two-level transition frequencies, and where 
the cavity can hence be adiabatically eliminated. In Sec. 
|V]we show that the squeezing, generated in the collective 
two- level degrees of freedom, can be released as a squeezed 
pulse of radiation, and we identify the field mode with 
the highest degree of squeezing. In Sec. VI of particular 
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relevance to solid state spin ensembles, we investigate 
the robustness of our squeezing and entanglement scheme 
towards inhomogeneities in the ensembles. Sec. VII 
concludes the paper. 



II. MODEL OF THE SYSTEM 

Assuming no inhomogeneity in the spin system, de- 
picted in Fig. la, and choosing h = 1, the Hamiltonian 
for the non-interacting cavity and spin degrees of freedom 
is given by 



A',, 



H Q = uJc+\f^w a 4> a + \Y^u b ai\ (1) 



where ui c is the angular frequency of the cavity mode with 
the annihilation (creation) operator c(cT), and the Pauli 




FIG. 2. A schematic picture of the system in which we im- 
plement spin squeezing. We have three harmonics oscillators 
of which one is inverted. The cavity oscillator operates as a 
bus for the transfer of excitations between the two ensemble 
oscillators. 



operators, a J z ' n , represent the j'th spin in ensemble n. The 
cavity mode magnetic field interacts with the spin degree 
of freedom, and with the rotating wave approximation, 
the interaction Hamiltonian reads 



A„ 



N„ 



(2) 



In the remaining text we assume equal interaction 
strengths gj >a = gjj = g\ for all j, and we will apply 
the Holstein-Primakoff approximation [2"4"] to describe the 
collective spin degrees of freedom as harmonic oscillators. 
The 6-ensemble is prepared with all spins in the spin 
down ground state, and we define the collective lowering 
operator b for the ladder of permutation symmetric states 
of the ensemble with 1, 2, ... flipped spins, 
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(3) 



As long as the ensemble has only few excited spins, the 
ladder of states is well described by the oscillator approx- 
imation, [6, w] ~ 1, and the dynamics under the ensuing 
interaction with the cavity field is readily solvable. This 
reveals the collectively enhanced coupling due to the large 
number Nj, of spins in the ensemble. 

The Holstein-Primakoff approximation applies for small 
permutation symmetric deviations from any collectively 
populated state |25j, and for the inverted a-ensemble, with 
all the spin prepared in the spins up excited state, we 
define the oscillator lowering (annihilation) operator 
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(4) 



The oscillator excitation number then counts the (small) 
number of spins flipped towards the spin down state, 
accompanied by the creation of photons in the cavity 
mode. 
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In terms of the oscillator ladder operators, our Hamil- 
tonian of the uncoupled systems writes 

Hq = u) c c'c — ujaO^a + LObtfb (5) 

and the interaction Hamilton is 

Hi = g a (a t c t + ca) + ffb (& t c + c t fe) (6) 

with g a = \/N a gi and g b = s/Nbgi, revealing the collec- 
tively enhanced coupling due to the large numbers, N a , 
and iVft, of spins in the ensembles. This Hamiltonian can 
be pictured as shown in Fig. [2j 

In the case of trapped three-level atoms or ions in an 
optical cavity, Fig. lb, we assume an off-resonant classi- 
cal laser field driving of the transition between the upper 
atomic ground state, labelled "spin up", and the opti- 
cally excited state, which in turn is coupled to the lower, 
"spin down" ground state via the emission of a cavity 
photon. With a sufficiently large detuning A, we may 
eliminate the optically excited state and thus retrieve a 
two-photon Raman process between the spin up and down 
states accompanied by the emission or absorption of cav- 
ity photons with a coupling strength g\ = fix/A, where 
f2 is the classical Rabi frequency, and x is the coupling 
strength (per photon) to the quantum field. Ensemble 
a, prepared initially in the spin up state, can undergo 
spin flip transitions and create cavity photons, while spin 
flips in ensemble b, prepared in the spin down state are 
accompanied by absorption of a cavity photon, and this 
system is thus also described by equation ^ and pj. 

The two set-ups in Fig. 1 involve different physical 
systems and fields at very different frequencies, and they 
offer different means of experimental control. The spin 
ensembles and the microwave cavity may be tuned in and 
out of resonance with each other thus turning their effec- 
tive coupling on and off, while the atomic ensembles are 
driven by a laser field, for which both the intensity and 
frequency offer means to control the coupling. Cavities 
for microwave and optical fields typically have different 
damping times, and while both spins and ground state 
atoms may have very long lifetimes, solid state spins may 
experience significant inhomogeneous broadening. We will 
include cavity damping throughout the following calcula- 
tions, and we will return to the effects of inhomogeneous 
broadening in Sec. VII. 

The Holstein-Primakoff approximation that we use is 
only valid as long as the excitation number in each en- 
semble is small. This will of course have to be checked 
in the calculations. The magnitude of squeezing can be 
qualitatively related to the number of quanta involved 
by recalling the dimensionless position momentum un- 
certainty relation (x 2 )(p 2 ) > 1/4 and the energy relation 
(p 2 /2 + x 2 /2) = (n) + 5, which suggests an amplitude 
squeezing factor ~ 1/ (n) . Spin squeezing by a factor 
10-50 thus involves flipping of a similar number of spins, 
which is indeed much less than the number of spins that 
we have in mind for the ensembles. 

The interaction Hamiltonian includes a term resembling 
the usual non-degenerate OPO Hamiltonian [2], and in 



a) b) c) 




20 40 60 20 40 60 20 40 60 



Nit Kjt tvt 

FIG. 3. The average number of excitations in the two ensem- 
bles calculated numerically with the parameters A = uj' — uj = 
10k, gb — 5k, and with g a = 0.6c?b , g a = 0.9g b , and g a = 1.2g b , 
in panels a, b and c, respectively. 

the absence of the 6-ensemble, we would expect the rather 
straightforward formation of a two-mode squeezed state 
of the cavity field and the a-ensemble. Since the cavity 
field leaks, this state will not live for long, and we hence 
couple the field to the second ^-ensemble and investigate 
if a long lived squeezed state can be created between the 
two ensembles. 

We first assume that g a ,9b > K, where k is the field 
decay rate. In this case the dynamics between the oscilla- 
tors is much faster than the damping of the cavity and 
for short times we can therefore neglect the decay term. 
For a simple and qualitative argument lets consider the 
resonant case, oj c = uj a = w&, where the time evolution 
operator in the interaction picture is given by 

U{dt) = e - lHldt « I - iHidt - ]^Hjdt 2 . (7) 

Applying ¥j\ to the initial state, we get 

W(dt)|0,0,0) = |0, 0, 0) - ig a dt\l,Q,l) 
dt 2 

- — (2^2,0,2) +^0,0,0) + ffaffb |l, 1,0)), (8) 

and we observe the last component with a simultane- 
ous excitation in both ensemble oscillators, characteristic 
for two-mode squeezing. The variance of the two-mode 
quadrature operator 

X ab {9) = \ (^(a + b) + ^( t + &t)) , ( 9 ) 

can be readily obtained, and for 9 = we get 

1 dt 2 

((AX ab (0)) 2 ) = - + —(3g 2 a Ag a g b ). (10) 

We clearly see a small reduction, provided g b > 3g a /4, 
where the asymmetry in the requirements on the coupling 
strengths is due to the respective OPO and beam-splitter 
like couplings to the ensembles. 
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FIG. 4. Variances of the collective quadrature observables 
X a b(0) ( re< 4 solid line), and X a b(n/2) (blue dashed line) for 
the coupling parameters and detuning, gb = k , A = 5k and 
g a = 0.5g b . 
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FIG. 5. The variance of X a b(0) for different detuning of the 
cavity. We have gt — 5k and g a — 0.9gt ■ The black (solid) 
line is for A = 75k. The red (dashed) line is for A = 50k. The 
blue (dotted) line is A = 20k and the green (dashed dotted) 
line is A = 5k. 



III. HEISENBERG-L ANGEVIN EQUATIONS, 
DAMPING AND STEADY STATE 



The bilinear Hamiltonian leads to simple Heisenberg 
equations of motion for the oscillator ladder operators. 
Cavity damping is modelled as a similar coupling to the 
continuum of free-space vacuum modes at one of the cavity 
mirrors, leading to a damping rate ft of the intra-cavity 
field accompanied by a white-noise term [2] , 



da 
dt 

~dt 
db 
~dt 

dtf 

~dt 
dc 

dt 
lit 



iua — ig a c* 
—iua' + ig a c 
—icob — ig b c 
iujtf + igbc' 



(k + iui')c — ig a a' — igbb - 
-(k — iui 1 )^ + ig a a + igbb' 



(11a) 
(lib) 
(11c) 
(lid) 
2ft Cm (<) (lie) 
2nc\ n {t), (llf) 



with [ci n {t), c| n (i')] = S(t — t'). In these equations we 
assume that the spin ensembles have degenerate excitation 
frequencies lu and the cavity field may have a different 
frequency u' . Starting from the oscillator ground states, 
all amplitudes have zero mean for all times, and the 
state is fully characterized by the second moments of the 
quadrature operators, which solve a linear coupled set 
of equations, following in a straightforward manner from 
(111. These equations are readily solved numerically, and 
in Fig. [3]we show the average excitation number of the two 
ensemble oscillators for different parameters. The mean 
excitation of the ensembles show two different regimes: 
One in which the excitation number undergoes a damped 
oscillation and one in which it increases exponentially. 

We can understand this behaviour if we define two new 
superposition mode annihilation operators 



d = (gba + g a b^)l \Jgl- gl 
e = (.gaa 1 + 9bb)/ijg%-gl, 



assuming gb > g a . It follows from (11) that 
d 



dt 



d = iujd, 



(12) 
(13) 

(14) 



and this "dark" mode d is uncoupled from the dynamics 
with a trivial time evolution. If we go to a rotating frame 
at the ensemble frequency uj, ^d — 0, then we can easily 
find the eigenvalues of the matrix that represents the 
homogeneous part of equation (111, 



Ai = 
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A 2 = 









(15b) 


A 3 = 
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2 ~ 




(15c) 


A 4 = 


i.A 
~2 ~ 


ft 

2 + 




(15d) 


A 5 = 


iA 


k 
~ 2 




(15e) 


A 6 = 


iA 
2 


k 
~ 2 




(15f) 



with G 2 



2iftA, where A 



w —ui. The homogeneous solutions can be written as linear 
combinations of e Xit , and for A = and gl > g^+ ft 2 , the 
real part of V G 2 is smaller than ft, ensuring a damped 
oscillatory behaviour of the solutions. 

In the parameter regime of a damped oscillatory so- 
lution, our system thus reaches a steady state, in which 
we will now investigate the squeezing properties. In Fig. 
[4] we show the variance of X a b(0) and X Q b(f), and we 
see that in steady state we obtain a squeezed state. The 
amount of squeezing depends on the detuning, the cou- 
pling strength and the decay rate of the cavity. In Fig. [5] 
we show the variance of X a b(0) for different values of the 
detuning and we observe squeezing for both short and 
large times. For large cavity detuning the damping due 
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FIG. 6. The minimally achieved variance for A = and 
different values of k and g a . 



The reduced dimensionality of the problem now allows a 
straightforward solution, and, e.g., for a(t) we have 



a(t) 



2 t 8 °" a i» 



9a - 9b 




A a{0) + w 



few- 



9b9b 



9l-9 2 b 



6t(0) 



9l-gl 



-^4.(0 



3b- 



(*-*') s 



5a56 f 



9a - 5b 



4(0 K- 



(19) 



Similar expressions can be found for , 6 and 6^ , and we 
can thus compute the variance of the quadrature X ab (0) 
observable, 



to cavity loss is small, and the spin ensembles show an os- 
cillatory dynamic between a significantly squeezed and a 
non-squeezed state, while for smaller detuning the damp- 
ing is more pronounced, and we quickly reach a squeezed 
steady state. In Fig. [6] we show the minimally achieved 
variance for different values of the coupling strengths 
to the inverted ensemble and cavity decay rate. Small 
damping and strong coupling lead to the largest degree 
of collective squeezing. 



((AX ab (0)) 2 ) = 
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9l + 9b 



iKt- 



8 y (5a + 9b) 2 

gl-: 

. 2g a g b (e t *+* 



+ e 



(.9a + 9bf 



/(*), 



(20) 



where f(t) is the contribution from the input field noise- 
terms. This is a double integral over the noise operators 
Ci n and c\ n , which by the commutator relations reduces 
to a single integral, which evaluates to 



IV. LARGE DETUNING - ADIABATIC 
ELIMINATION OF CAVITY FIELD 



If the cavity field detuning is large, A = uj' — oj g a ,gb, 
we may eliminate the field operators from the problem by 
setting the derivative to zero 

d c 

= — = -(« + iA)c - ig a a) - ig b b - v2Kc in (t) (16) 
at 

with the solution 



c(t) 



.g a a^ +g b b 
■i r~7 — 

K + lA 



'2k 



iA 



(17) 



This expression can be inserted into the Heisenberg- 
Langevin equations for a and b such that we get 



da 
~di 
da) 
lit 
db 
~dt 
d&t 

~dt 



9 2 a , 

^r a ^ 

k — iA 

K + lA 



9*9b ^ 

K — iA 

g a gb , 



ig a y2K 

k — iA 
ig a V2~K 



k + iA K + iA 



Cin (^) 



(18a) 
(18b) 



9?, 



, g a gb + igbVZn ,,, , 1Q , 

a-0 i~ a ir c in\t) (18c) 

K + iA K + iA K + iA v ; 



._L 6 t 
k — iA 



g a gb . igbVZn t 



iA 



iA 



/(*) = 



2k/ 



- 1 (5a - 56 ) 



8 (5a + 5&) 



(21) 



In Fig. [7] we have plotted equation ( 20 1 for different cou- 
pling strength and wee see that the behaviour resembles 
the result from the numeric analysis. 



The exponential terms in equation ( 20 ) and in the 



integral over the noise terms are decaying with time, if 
56 > 5a 7 an d in the long time limit the variance converges 
to 



((AX ab (0)f 



9l 



%a+56) 2 ' 



(22) 



which also confirms the squeezing results. 



Furthermore we know that the exponents in ( 20 ) are 



small in the adiabatic regime, which enables an expan- 
sion of the exponential functions such that we get an 
approximate expression for the variance, 



((AX ab (0)) 2 



1 



9 ? 
9 a- 9 b 



Kt 



(23) 



which is valid for small times only. 

If we choose a finite interaction time t, such that 



t' 



A 



±7r, and if A > k, the complex exponential 



terms in ( 20 ) undergo a complete change of sign before 
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FIG. 7. The collective quadrature variance, determined by 
adiabatic elimination of the detuned cavity field, (20 1, with 
A = 75k and gt — 5k. We vary the coupling strength, while 
remaining in the adiabatic regime. For the black (solid) line 
we choose g a = 0.9gb, corresponding to the black line in 
Fig. pM For the red (dashed) line we choose g a = 0-Sgb, the 
blue (dotted) line g a = 0.7gr& and the green (dashed dotted) 
g a = 0.6g b . 



any appreciable damping of the exponential functions 
occurs. To first order in # this leads to a minimum in 



the variance 



((AX ab (0))' 



(ga ~ 9b) 2 

% a + g b )'- 



+ 



2 5a(56 -9a) K 



8(g a +g b ) 2 A 



(24) 



which approaches zero for g a approaching g b and ^ ap- 
proaching zero. These minima are readily observed in 
Fig. 



annihilation operator cx g a a + g b b couples directly to the 
cavity field. Ideally, the squeezing properties in this mode 
should be converted into cavity photons leaving the cavity 
at a rate n and thus forming a finite pulse with favourable 
squeezing properties. 

The Hamiltonian coupling and the damping and input 
noise fields leads to Heisenberg-Langevin equations for 
cavity field and the coupled mode, that we can write on 
the form 



-r(t) = M in r(t) 
at 



where 



and 



r(t) 



g a a(t) + g b b(t) 
c(t) 



r in (t) 



Cin{t) 



(27) 



(28) 



(29) 



and the matrix M in contains the coupling terms from the 
Heisenberg-Langevin equations. 

The input/output formalism 0|SB], relates the input 
and output fields by the equation 

c ut(t) - c m (t) = V2nc(t) 

This allows us to write the Heisenberg-Langevin equa- 
tions in an alternative form, driven by the output noise 
operators, 



V. SQUEEZING IN THE OUTPUT FIELD 

Conventional squeezed light sources are based on op- 
tically non-linear materials and they are subject to a 
competition between the non-linear processes and losses. 
In our scheme, the squeezing of the spin ensembles is 
not due to a non-linear interaction but to an exchange of 
quanta with an initially inverted medium, and our analy- 
sis suggests that a potentially high degree of squeezing is 
obtainable. This squeezing may be relevant for entangle- 
ment operations on the spin degrees of freedom, but they 
may also be coupled resonantly to the cavity field, and 
thus lead to production of a squeezed output field from 
the cavity. To accomplish this mapping, we will apply a 
7r-pulse to the inverted ensemble, such that its excitations, 
correlated with the ones in the other ensemble, can be 
converted into photons according to the Hamiltonian of 
the system, where the annihilation and creation operators 
of ensemble a now attain the same roles as in ensemble b, 



H' a = uj c c} c + u) a a) a + ui b b^b 

Hi = g a (a f c + c)a) + g b (tf c + c f 6). 



(25) 
(26) 



We observe that the spin ensembles also here display an 
uncoupled collective mode, while the collective mode with 



dt 



r(t) = M out r(t) - V2nr out (t), 



where 



?out(t) 





Cout{t) I ' 



(30) 



(31) 



and M out differs from M in by a change of sign of the 
damping term of c. 

The Laplace-transform turns the differential equations 
into algebraic forms 



sf (s) - r(0) = M m r(s) - V2nr m (s) (32) 
= M out r(s) - V2~Kr out (s), (33) 

with t = at the time of the 7r-pulse. From this we can 
eliminate f (s) and obtain the expression 



?cmt(s) =(s - M out )(s - M in ) 

+ r(0l 
v/2kT 



_! (V2KT in (s) - r(0)) 
V2k 



(34) 



By carrying out the matrix inversion and mulitplications, 
and performing the inverse Laplace-transform we finally 
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obtain the expression 

a(t)(g a a(0) + g b b(0)) - /3(t)c(0) 



o.2o r 



2 k 



+ / P{t-t')c in {t')dt' + c in (t), (35) 
Jo 



with the functions 

w = ^! e -lc«+iA) isinh (i m ) 



(36) 



/3(t) = 2 K ( / ^-^sinh(-m) -cosh( 1 m))e-5(-+^ 

(37) 



where = \J k 2 



\g\ - A 2 - 2ikA. With the ex- 



pression for the output-field in terms of operators with 
known noise properties, we can compute the degree of 
squeezing at any moment of time and non-classical corre- 
lations between different times. Since the field is emitted 
in a pulse after the inversion of the a ensemble, we expect 
that a time dependent mode function for the field u(t) 
exists, for which the squeezing is maximal. We define the 
corresponding mode operator [27] 



u(t)c out (t) dt 



where u(t) is normalized such that 

\c u , 41 = i. 



(38) 



(39) 



The noise in the corresponding field quadratures is mea- 
sured by homodyne detection with either a local oscillator 
with the appropriate time dependence, or it is determined 
by a weighted temporal integral of the homodyne signal 
obtained with a constant local oscillator. The degree of 
squeezing of course depends on the shape of the chosen 
mode function u(t), and we have studied different candi- 
date mode functions, including a simple exponential decay 
with rate k, the real square root of the mean intra cavity 
photon number .251. and a general complex function u(t) 
which is numerically optimized to yield the largest degree 
of squeezing. 

The variance of X Cu ( ~ ) for these three choices of mode 
functions are shown in Fig. [HJ and examples of their time 
dependence are illustrated in Fig. [9j We note that if 
the emitted light is squeezed within a single mode, the 
choice of a different mode function is equivalent to the 
admixing with a vacuum state with a standard variance 
of \. The exponentially decaying function thus has a 
too small overlap with the optimum mode to show any 
significant squeezing. For small damping, the square 
root of the emitted power provides a good ansatz for the 
squeezed light mode. For stronger cavity damping the 
numerically optimized mode deviates in particular for 
early times and shows significantly better squeezing. The 
dependence of the optimal squeezing on kappa in Fig. 
[8] is similar to the spin squeezing dependence shown in 
Fig. [6j and suggests that the squeezing properties are 
successfully transferred to the emitted light. 
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FIG. 8. The output field quadrature variance with u(t) = 
Ae' Kt/2 (blue dashed dotted), u{t) = B^c)) 1 ' 2 (red dashed) 
and a optimized u{t) (solid black line) for different values of 
k. We take A c = O.Olgt, and g a ~ 0.9gt- 




FIG. 9. Non-normalized u{t). We have u{t) = Ae~ Kt/2 (blue 
dashed dotted), u(t) = A({c^ c)) 1 ^ 2 (red dashed) and optimized 
u(t) (solid black line). We take g b = 1.25k, A = 0.001k and 
g a = 0.9g b . 



VI. SQUEEZING WITH INHOMOGENEITY 

For solid state spin ensembles, it is not a good ap- 
proximation to assume perfect degeneracy of the spin 
excitation energy. Depending on the spin system itself 
and on the purity of the host material, individual spins 
have excitation frequencies that may vary within a MHz 
wide profile or more, and this both detunes the spins from 
the cavity and causes precession of the spin excitation at 
different frequencies, leading effectively to a damping of 
the quadrature correlations. 

We model the inhomogeneous broadening in the spin 
ensembles by allowing the energies to be different for each 
spin, 



H = ujc + u ha al' a + \Y, u 3^ b (40) 



and they may also experience different interaction 
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FIG. 10. The variance of X o (,(0) for different broadening 
width. Defining gb = \/^2 k |<?fc,&| 2 , we have at the top <?t, = 5k, 
= , a. | 2 = 0.9<76 and A c = 75k and in the bottom 

we have the same but A c = 0.5k. The black (solid) line is for 
a width of Ogb- The red (dashed) line is for 0.05k. The blue 
(dotted) line is 0.5k and the green (dashed dotted) line is for 
5k. 



strengths with the quantized field mode 

N a N a 



Hi = 9jA a + ac + a - a ^) + X! 9j,b{ <J +' c + a - 



(41) 



To use the Holstein-Primakoff formalism we group the 
spins in frequency intervals with nearly identical frequen- 
cies, for which the harmonic oscillator operator represen- 
tation may be used for the a-ensemble, 



1 



(42) 



and for the 6-ensemble, 




(43) 



With the inverted a spins, we now have the uncoupled 
and the interaction Hamiltonian 



Hq = UcC^c — ujjOjdj + cjjbjbj (44) 



Passing to a rotating frame, and introducing detunings 
Aj and A c with respect to a central spin frequency, we 
obtain the Heisenberg-Langevin equations 



dcij . t 

— = lAjdj - ig j>a & 



dt 
dh 
dt 

d A 

dt 
d c 



= -iAjdj + igj, a t 
— iAjbj ~ igj t bC 



—iAjbj + igjfiC 

— (k + iA c )c 

~ //./■■'"■ + 9j,bbj) - V2Kc in (t) 

3 

— {k - iA c ) 



(46a) 

(46b) 
(46c) 

(46d) 
(46e) 

(46f) 



This leads to deterministic linear equations for the sec- 
ond moments of the quadrature operators, and in Fig. 
[10] we show results for the variance of the collective spin 
mode with different values of the inhomogeneous broad- 
ening. In the calculations we have assumed a Gaussian 
distribution of the spin frequencies. The upper part of 
the figure shows the case where the cavity is detuned with 
respect to the spin frequency distribution - the squeezing 
occurs on a long time scale, and the precession of spins is 
likely to be the cause of the significant reduction of the 
squeezing already for moderate broadening. The lower 
part of the figure shows the almost resonant coupling to 
the cavity field and a resulting faster squeezing, which is 
also observed for weak inhomogeneities. In addition, the 
cavity coupling itself may effectively separate the collec- 
tively coupled spin mode from the other modes and thus 
suppress the effect of inhomogeneities [29\ 130], 



VII. CONCLUSION 

We have shown that we can create two-mode spin 
squeezed and entangled states of separate spin ensem- 
bles in a microwave cavity and atomic ensembles in an 
optical cavity from a state where one ensemble is inverted 
and excitations are transferred between the ensemble via 
the cavity field. We have investigated this process with 
Heisenberg-Langevin equations, which readily provide so- 
lutions for the state of the spins and the cavity field, and 
in the case of a large cavity detuning, we have provided 
simple expressions for both the long time limit (22 1, and 



for optimal interaction times where maximum squeezing 
is obtained ( 24 ) . These results could be investigated in 



existing systems with trapped atoms and ions in optical 



9 



cavities, and we have extended the theory to include inho- 
mogeneity in the ensembles, which may be of relevance for 
electron spin ensembles in solids. Inhomogeneities reduce 
the squeezing for long times and change the optimum 
values, but should still be observable, e.g., in experiments 
with NV centers, and possibly even with species of differ- 
ent particles. We have also shown how the intra cavity 
squeezing can be transformed into a squeezed pulse of 
radiation, emitted on demand by the system, and be- 
ing thus of interest for quantum communication schemes, 



and possibly also for studies of the interactions between 
quantum radiation and superconducting qubits 31j. 
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